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We perform a combined analysis of B — > tttt and B — > irK decays with the current ex- 
Qh. perimental data. Assuming SU (3) flavor symmetry and no new physics contributions 

*^ to the topological amplitudes, we demonstrate that the conventional parametrization 

in the Standard Model (SM) does not describe the data very well, in contrast with 

X 

a similar analysis based on the earlier data. It is also shown that the introduction 
of smaller amplitudes and reasonable SU (3) breaking parameters does not improve 
the fits much. Interpreting these puzzling behaviors in the SM as a new physics 
(NP) signal, we study various NP scenarios. We find that when a single NP am- 
plitude dominates, the NP in the electroweak penguin sector is the most favorable. 
However, other NP solutions, such as NP residing in the QCD-penguin sector and 
color-suppressed electroweak penguin sector simultaneously, can also solve the puzzle. 
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1 Introduction 



In the Standard Model (SM), rare non-leptonic decays B — > tttt and B — > ttK provide 
valuable information on the inner angles of the unitarity triangle of Cabbibo-Kobayashi- 
Maskawa (CKM) matrix, and have been widely studied. For this purpose the measurement 
of time-dependent CP-asymmetry given by 
r(g(t) -» fop) - T(B(t) -» /cp) = 

r(B(t) - /cp) + r(s(f) - /cp) v ' v ' w 

is essential. Here Acp and Scp represent direct and indirect CP asymmetries, respectively. 

Specifically, B — > 7T7t decays measure the angle a through the isospin analysis [1]. 
The information on 7 can be obtained from B — > 7T.K" data [2, 3, 4]. In addition, if 
a new physics (NP) beyond the SM exists, it can significantly affect these processes by 
contributing to penguin amplitudes. Therefore these decay processes are also a sensitive 
probe of NP [5, 6, 7]. 

Assuming i) SU (3) flavor symmetry of strong interactions and ii) smallness of annihi- 
lation and exchange topologies, Buras, Fleischer, Recksiegel and Schwab [5] concluded the 
B — > tttt and B — > ttK data strongly suggest a NP in the electroweak penguin sector of 
B — > ttK decay amplitudes. On the other hand, Chiang, Gronau, Rosner and Suprun [8] 
demonstrated that the x 2_ fit to the same data does not show any significant deviation from 
the SM. It should be noted that Buras, et.al. [5] assumed that there is no significant NP 
contribution to B — > tttt decays, and used some of B — > tttt data with small experimen- 
tal errors to predict the hadronic parameters of B — > ttK amplitudes. Chiang, et.al [8] 
included all the available B — > tttt and B — > 7rfC data into their fit. 

We considered B — > 7rfT data only [6] , and showed that the SM fit faces some difficulties 
and NP in the electroweak penguin sector is strongly favored in accord with [5]. We 
concluded that the discrepancy between [5, 6] and [8] is due to the dilution of NP effects 
by including all the data in [8]. 

In this paper, we perform x 2 fitting to the current data in the SM and also in the 
presence of NP. We show that even if following the approach of Chiang, et.al [8], i.e. the 
X 2 fitting with all the available data, we get much worse x 2 fit than in [8]. We calculated 



Ax 2 -the contribution of each data point to x 2 value- to trace the source of this puzzling 
behavior. To improve the fit in the SM, we introduce i) smaller amplitudes, and/or ii) 
reasonable S'L r (3)-breaking effects to the fits. It turns out that these corrections do not 
solve the puzzle satisfactorily. 

Interpreting these difficulties in the SM fits as a NP signal, we introduce NP parameters, 
such as, a new weak phase in the amplitudes. We consider various NP scenarios. We 
introduce three types of NP, NP in the electroweak penguin, color-suppressed electroweak 
penguin and QCD penguin. When a single NP amplitude dominates, NP in the electroweak 
penguin is the most favorable solution, supporting the findings in [5]. A given specific NP 
model, however, contributes to all the NP amplitudes in general. In light of this we also 
considered the possibility two or more NP amplitudes are enhanced simultaneously. 

The paper is organized as follows. In Section 2 the SM fittings are considered. In 
Section 3 we perform various NP fittings. The conclusions and discussions are given in 
Section 4. 



2 SM fitting 

The topological amplitudes provide a parametrization for non-leptonic 5-meson decay pro- 
cesses which is independent of theoretical models for the calculation of hadronic matrix 
elements [9]. The decay amplitudes of B — > 7T7r's which are b — > dqq (q = u, d) transitions 
at quark-level can be written as 

y/2A(B+ -> tt+tt ) = -(T + C + P EW + Pg w ), 

A(B°^n + n-) = - (t + P + |p£ w + E + PA^j , 

V2A(B°^n°n°) = - (c - P + P EW + ^P° w - E - PA^. (2) 

Here T, C, P, 

^EW' & an d PA represent tree, color-suppressed tree, QCD-penguin, (color- 
suppressed) electroweak-penguin, exchange and penguin annihilation diagrams, respec- 
tively. Similarly, B — > ttK decays which are b — > sqq (q = u, d) transitions at quark-level 
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are described by 
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where primes indicate b — > s transition. The corresponding decay amplitudes for the 
CP-conjugate modes can be obtained by changing the sign of weak phases while keeping 
CP-conserving strong phases unchanged. 

We can further decompose the QCD penguin diagrams, P and P', depending on the 
quarks running inside the loop, 



P = V ud V* b P u + V cd V* b P c + V td V* b P t 

= V ud V* b (P u - P c ) + V td V t t(P t - P c ) 

= Puc&^ ~\~ Ptc& ^ 1 

P' = V us V* b P' u + V cs V* h P' c + V ts V* b P' t 

= v us v: b {p' u - p') + v ts v*(p; - p' c ) 

= K^-pL (4) 



where we have used the unitarity relation for CKM matrix elements and explicitly written 
the weak phase dependence for the amplitudes. These notations and conventions will be 
used throughout the paper. We can estimate the relative sizes of the amplitudes based on 
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the color-, CKM-, and loop-factors, 
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where A is expected to be order of 0.2 ~ 0.3. We will call the decay amplitudes parame- 
terized as in (2) and (3) and the hierarchy in (5) the conventional parametrization in the 
SM. 

The decay amplitudes containing only dominant terms, T^'\ 2 ,PfJ and P^ w are 
given by 3 

V2A(B+ -> tt+tt ) = -(T + 0)e% 

A(B°^tt + tt-) = -(Te^ + Pe^), 
V2A(B° -> ttV) = -O l7 + Pe" i/3 , 

A(P + ^7r+K°) = -P', 

V^A(S + -> 7T ^+) = P' - TV 7 - OV 7 - P^ w , 

A(P°^7T^ + ) = P'-TV 7 , 

v^P -> 7T°K ) = -P' - OV 7 - P^ w . (6) 

Here, we have written P^ as P^ for the simplicity of notations. 

The current experimental data for the CP-averaged branching ratio (BR), the direct 
CP-asymmetry (Aqp) and the indirect CP-asymmetry (Sep) are shown in Table 1. It 
immediately shows some puzzling behaviors which are difficult to understand if we believe 

2 We include C to the amplitudes, although according to (5) it is subdominant. Otherwise, we get 
extremely poor fit [6]. We hope that this problem will be solved within the SM framework. 

3 We may also think of T(') and CO as + pQ and - pQ , respectively. See [8, 6], for details. 
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1V1U Qc 


-DiXllU ; 


A t — 1 T — > 

^*CP 


•S'cp 


B + -> 71 + 71° 


5.5 ±0.6 


0.01 ±0.06 




B° -> 7T + 7T- 


5.0 ±0.4 


0.37 ±0.10 


-0.50 ±0.12 


£?° -> 7T 7T 


1.45 ±0.29 


0.28 ±0.40 




5 + -> 7T + K° 


24.1 ± 1.3 


-0.02 ±0.04 




£? + -> 7T°K + 


12.1 ±0.8 


0.04 ±0.04 





B°^n-K + 18.9 ±0.7 -0.115 ±0.018 

B°^n°K 11.5 ± 1.0 0.02 ±0.13 0.31 ± 0.26 

Table 1: The current experimental data for CP averaged branching ratios (BR), direct CP 
asymmetries (A CP ) and indirect CP asymmetries (S'cp) for B — > rnr and B — > ttK decays 
filed by HFAG [12]. 

(5) and (6). Firstly, (5) suggests that BR(B° — > tt°-k ) should be about 3 times lower than 

the data (B — > tttt puzzle) [7, 13]. Secondly, the ratios 

2BR(B + -> n°K + ) 
c = BR{B+ -> ir+K°) ' 
-» vr-^+) 

n ~ 2BR(B° ^n°K°) { ) 

should equal to a good approximation. However, the data shows about 1.5a difference. 
We should say that this so-called R c jR n problem is not so statistically significant now. 
Thirdly, we expect from (6) that 

A CP (B + -> 7r°K + ) « A CP (B° -> 7T-X+). (8) 

The data deviate from this relation by about 2.7er level. Finally, the dominant terms in 
Scp(B° — > ti° gives sin 2/3 which is quite precisely measured from b — > sec modes to 
be sin 2/5 = 0.685 ± 0.032 [12]. The current data shows about 1.43a difference. These last 
three are usually called " B — > irK puzzle" [14]. 

Assuming the exact SU (3) flavor symmetry, we can relate the topological amplitudes 
of B — > 7T7T decays to the corresponding amplitudes of B — > 7T.K" decays as follows: 

T £ = 
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\C'\ (eV) 

| W | V ) 


3 82 ± 48 


3.20 ± 0.52 


4.02 ± 0.44 


5 C ' 


-40.2° ±9.7° 


-14.3° ± 18.0° 


-47.0° ± 10.5° 


\P'\ (eV) 


48.9 ±0.7 


47.2 ± 1.7 


36.9 ±4.8 


l^cl (eV) 






24.1 ±6.4 


SpL (eV) 






178° ± 2° 



Table 2: Results for "SM fit I", "SM fit II" and "SM fit III". See the text for details 
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In addition, it is known that the Wilson coefficients for the electroweak penguins c-j and 
Cs are much smaller than eg and cio [10] in the SM, which leads to a relation between the 
electroweak penguin diagrams and trees in the S77(3)-limit [11], 



p ew = 3 /-^R(T' + c') + 3 - c -^^R(r-c>), 

4 d ± c 2 4 ci - c 2 

3 c 9 + c 10r{t , + _ 3^-^o _ 



EW 



4 ci + c 2 4 ci - c 2 

Here, R is given by a combination of CKM matrix elements, 
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Using these SU(3) relations, we have 6 parameters to fit in (6): |T'|, |C"|, \P'\, two 
relative strong phases and 7 (SM fit I). Now we can perform the fit to the current exper- 
imental data which are given in Table 1. Since in -P^w 1S neglected in (6), we use for this 
fit 



-Pew - - 9 10 R(T' + C), P E ^ - 0. 



2 ci + c 2 



(12) 
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The inner angle (3 of the unitarity triangle is strongly constrained and is given by sin 2(3 = 
0.725 ± 0.018, so we fixed (3 = 23.22°. Other parameters used as inputs for the fit are as 
follows: A = 0.226, c x = 1.081, c 2 = -0.190, c 9 = -1.276a cm , c w = 0.288a cm . The result 
for this fit is shown in Table 2. 4 It can be seen that the SM with exact SU(3) symmetry 
gives Xmm/dof = 18.8/10(4.3%), which is quite a poor fit. To trace the observables which 
make the fit poor, we list Ax^ lin -the contribution of each data point to the Xmm m Table 3. 
From Table 3 we can see that the observables, BR(B° -> ir°K°), A C p(B+ -> n°K+) and 
Scp(B° — > ir°K°) which caused the B — > nK puzzles are exactly those with large A^mm- 
They are about 1.8a ~ 2.2a away from the best fit values. 

An alternative way to see the discrepancy between the SM and the experiments is 
to remove the observables which give large Ax^ in from the fit and predict them from 
the fitted parameters of the remaining observables. For example, we dropped the data 
for BR(B+ -> n+K°), BR(B° -> n~K + ), BR(B° -> 7i°K°), A CP (B+ -> tt°K+) and 
S CP (B° -> n°K°) whose A^^'s from "SM fit I" are greater than 1 from the x 2 fitting. 
The results for this approach (SM fit II) are shown in the 2nd columns of Tables 2 
and 3. From Table 2, we can see the quality of fitting has improved dramatically while 
the values of parameters are consistent with those of "SM fit I" . Also Table 3 shows that 
all the observables considered are excellently described by the SM parametrization. Now 
we can predict the omitted observables from the fitted values in Table 2. The predictions 
(deviation from the best fit values) for BR(B + -> n + K°), BR(B° -> n-K + ), BR(B° -> 
ir°K°), A CP (B+ -> tt°K + ) and S C p(B° -> n K°) are 21.0 ± 0.6 (2.1a), 18.5 ± 0.4 (0.46a), 
8.17 ± 0.16 (3.3a), -0.065 ± 0.002 (2.6a) and 0.81 ± 0.0001 (1.9a), respectively. These 
deviations imply that the B — > nK puzzles are more serious than the estimations given 
below (7) and below (8). 

Until now we have assumed exact SU(3) flavor symmetry. Before considering NP as 

a solution of these B — > nK puzzles, we proceed to improve the SM parametrization by 

4 We have also checked the case where P' EW ~ — § ^^+02) ^ EW ^ ^" anc ^ ^ EW ^ s ^^ted as in [8] . In this 
case we obtained <?ew = 0.36 ± 0.33 which is far away from the SM expectation <5ew = 1- Therefore we 
used the exact SU(3) relation (10) for this fit. 
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Observable SM fit I SM fit II SM fit III 



7~> 7"~) { 7~)-l- 
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1.2 




0.17 
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-> 7T°^ + ) 


0.025 
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0.59 
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1.3 




0.98 


BR(B° - 




4.8 




1.4 



Acp(B+ - 


-> 7T + 7T ) 
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0.028 


0.028 


AcAB - 
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9.9 x 10~ 5 


0.10 


1.1 


A CP (B° - 


-> 7T 7T ) 


0.50 


0.029 


0.40 


A CP (B+ - 


-> TT + iT ) 


0.25 


0.25 


0.18 


A CP (B+ - 


■* 7T ^+) 


3.1 




3.1 


A CP (B° - 


» 7T-^+) 


0.68 


0.044 


1.2 


A CP (B° - 


-> TT°K°) 


0.33 


0.15 


0.49 



S C p(B° - 7T+7T-) 


0.85 


0.00013 


0.19 


<? CP (P ^ 7r°X°) 


3.9 




5.9 



^ nin -the contribution of each data point to the Xmin- 



including smaller amplitudes we have neglected in the above analysis and/or by taking 
SU(3) breaking effects into account. First we include P' uc and P^w which are subdominant 
according to (5) in the decay amplitudes. Then we the decay amplitudes of B — > nn and 
B — > 7ri^ are corrected to be 

= -{T + C)e l \ 
-- - (re* + Pe~ if> ) , 
-- -C^ + Pe- i(} , 

-P' --P' C 4- P' 



V24(£+ - 
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v^A(S° - 


- ttV) 


A(B+ - 




V2A(B + - 


> 7T°^ + ) 



^ - ^ew " ^ew " (t' + C" + P; c ) e* 
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V2A(B° - n°K°) = -P'-P^- 1 -P^-[C'-P' U ^\ (13) 
We also incorporate the factorizable SU (3) breaking effect to the tree amplitude [8] so that 

where f n (fx) is the decay constant of n (K). For numerical analysis we used f n (K) = 
131(160) (MeV). For color-suppressed tree and QCD penguin amplitudes we still use the 
SU(3) relation (9). We also use the relation (10) for electroweak penguins in terms of trees 
(SM fit III). 

As can be seen in Table 2, the xLm does n °t improve at all. In addition |P^ C | which 
should be much smaller compared with \T'\ does not follow this hierarchy. We also see 
that Ax^'s for S C p(B° -> vr ^ ) and A CP {B + -> tt°K + ) in Table 3 are still troublesome. 
Therefore we conclude that the inclusion of factorizable SU (3) breaking effect in and 
smaller amplitudes P' uc and P^w a l° ne does not help improving the SM fit. 

Now we consider the effect of reasonable SU(3) breaking. To do this we introduce two 
parameters be and bp to represent the SU(3) breaking for the color-suppressed tree and 
QCD penguin so that 

C , V ud P V td 

c>- hc v:; p>- bp V t 

We added two free parameters b c and b P to "SM fit III" (SM fit IV) and obtained 
be = 3.5 ± 6.6 and bp = 1.7 ± 0.7 (Xmm/^ / = 4.7/5). Although they have huge errors, 
the central values require too large SU(3) breaking effect, considering the fact that it is 
expected to be at most 20 — 30%. To make matters worse, not only \P' UC \ is too large but 
7 = 41° ± 5° is much lower than that obtained in the global CKM fitting [15]. 



(15) 



3 NP fitting 

We have seen in Section 2 that the SM parametrization does not describe the experimental 
data very well. Although the discrepancy is about 2 — 3 a level and we cannot rule out the 
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NP fit I 


NP fit II 


NP fit III 


Xmm/dof (quality of fit) 


6.28/7 (51%) 


7.92/7 (34%) 


8.3/7 (31%) 


7 


71.7° ± 5.7° 


71.1° ±8.4° 


53.2° ±8.7° 


\T'\ (eV) 


5.21 ±0.27 


5.23 ±0.28 


5.40 ±0.30 


8t> 


30.3° ±5.5° 


32.9° ± 13.6° 


75.0° ±40.9° 


\C'\ (eV) 


3.25 ±0.54 


3.56 ±0.54 


4.41 ± 0.43 


5 C > 


-14.5° ± 18.7° 


-24.4° ± 15.4° 


-0.4° ±35.8° 


\P'\ (eV) 


48.6 ±0.7 


47.4 ±4.3 


25.2 ±8.0 


5np 


7.6° ±4.3° 


-5° ± 2° 


-100° ± 44° 


-Pew.npI ( e V) 
0ew (eV) 

-Pew.npI ( e ^) 
0ew (eV) 
I^npI (eV) 
0p (eV) 


20.1 ±4.7 
-87.4° ± 4.5° 


33.6 ±26.7 

-88° ± 3° 


39.0 ±19.3 
-1.94° ±2.12° 



Table 4: Results for "NP fit I", "NP fit II" and "NP fit III". See the text for details. 

SM yet, it would be interesting to investigate whether a new parametrization coming from 
NP will improve the fitting. 

Since the parameterizations in (6) and (13) can perfectly fit to the B — > nn data, we 
will assume that NP appears only in the B — > nK modes. A given NP model can generate 
many new terms in the decay amplitudes with their own weak phases and strong phases. 
To simplify the analysis we adopt a reasonable argument that the strong phases of NP are 
negligible [6]. With this assumption it is enough to introduce just one NP amplitude with 
effective weak phase for each topological amplitude. We assume there is no NP contribution 
to tree amplitude V and color-suppressed tree amplitude C Then the decay amplitudes 
can be written as 

V2A(B + -> tt+tt ) = -{T + C)e l \ 
A(B°^n + n-) = - (Te^ + Pe~^) , 
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Observable NP fit I NP fit II NP fit III 



BR(B^ - 


+ 0\ 
-> 7T + 7T U ) 


3.5 x 10 


0.071 


0.037 


BR(B° - 


_i_ \ 

-> 7T^7T j 


0.056 


0.086 


0.030 


BR(B° - 


-> 7T 7T ) 


2.4 x 10~ 7 


0.21 


0.032 


BR(B + - 


-> 7T + K°) 


1.7 


0.13 


0.089 


BR(B + - 


+ 7T°^+) 


0.16 


0.29 


0.70 


BR(B° - 


> TI-K+) 


1.1 


0.026 


0.97 


BR(B° - 


+ TT°K°) 


0.099 


0.46 


1.2 





-> 7T + 7T ) 


0.028 


0.028 


0.028 




— > 7T + 7T~) 


0.14 


0.31 


0.023 


A C p(P° 


-> 7r°7r°) 


0.021 


0.037 


0.48 


Acp(£+ 


-> 7T + K°) 


0.25 


2.0 


0.024 


A CP (P+ 


-> 7T°^+) 


0.13 


0.42 


0.15 


A CP (P° - 


-> 7T-^+) 


0.16 


0.29 


0.097 


A CP (P° 


-> 7T°^ ) 


1.8 


0.0013 


1.4 



S CP (P° - 7T+7T-) 


0.18 


0.058 


0.042 


S C p(P -7r°K ) 


0.49 


3.5 


2.9 



^ nin -the contribution of each data point to the Xmm- 



V2A(B° -> ttV) = -C* 7 + Pe~ i/3 , 
^,4(5+ -> 7r° J ft: + ) = P' - TV 7 - CV 7 - P EW - P^pe^ p - P E1 ' ■ u 



EW,NP C 



g-'EW.NP 6 J 

A(B°^n-K + ) = P'-TV 7 -P^pe^-^P E ^ )NP e^w, 
>/2A(B° - 7r°K°) = -P'-C'e^-P^ + PV^-^EWNP 1 



-^E^NPe^-. (16) 

Note that we included NP contribution to the color-suppressed electroweak diagram. This 
is because this contribution need not be suppressed compared to the electroweak penguin 
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while it is actually suppressed in the SM. 

This description of NP has 7 additional parameters, overall strong phase <5np relative 
to that of P' which we set to be zero, three real NP amplitudes and three NP weak phases. 
Using all the new parameters in fitting makes statistics quite poor. So, at first, we assume 
one NP terms dominate and neglect the others. 

First, we consider only the effect of P EW (NP ^ I) which corresponds to the solution 
considered in [5]. Table 4 shows that we obtained an excellent fit for this scenario. In 
addition, we can see in Table 5 that all the puzzling behaviors of the SM have disappeared. 
The largest deviation from the best fit parameters is at most 1.3a. 

Now we consider a scenario where P E w,np dominates (NP fit II). We can see from 
Table 4 that this fit is also acceptable. However, the data for Scp(B° — > ir°K°) is a little 
bit away from the best fitted values. 

Similarly we can consider the case where only P^ p exists (NP fit III). In this case, 
although Xmm/dof is acceptable, Axf nin of Scp(B° —> tt°K°) is not so satisfactory. 

We can analyze more general cases of having two-types of NP simultaneously. We have 
1 1 parameters to fit in these cases. For each case we obtained several acceptable solutions, 
which is due to the low statistics. With only P EW np an< ^ ^ew np (NP fit IV) we get two 
distinctive solutions in Table 6. It is interesting to note that the best solution of "NP fit IV" 
favors the enhancement of both P EWj np an d Pew.np contrary to Ref. [5]. As mentioned 
earlier, P E w,np ^ s n °t necessarily color suppressed and can be as large as P EWNP . The 
second solution corresponds to that found in Ref. [5], i.e. NP in the electroweak penguin 
sector. 

The scenario of having non- vanishing P^ P and P E w,np (NP fit V) shows that P^ P need 
not be suppressed and can be almost as large as P E w,np- Even when P E w,np vanishes (NP 
fit VI), the large contribution of P^p and Pewnp can solve the B — > nK puzzle. 

Now we consider the simultaneous contribution of all the possible NP contributions, i.e. 
Pnp,Pew,np and P E w,np f° r completeness (NP fit VII), although we have poor statistics 
for definite prediction. As expected, we obtained many physically acceptable local minima. 
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IN JT 111 IV 


NP fit V 

1NJT 111 V 


NP fit VT 

IN 1 111 VI 


XLJdof (quality of fit) 


2.45/5 (78.0%) 

9 79 (7A A°7r>\ 
Zi. t Zij o yl 4.4 /o J 


4.55/5 (47%) 


0.51/5 (99%) 
u.y i j o yi> i /oj 


7 


69.7° ±5.6° 

( 0.0 IE 0.0 


62.4° ± 8.4° 


60.7° ±8.8° 

OO. U x (.1 


\T'\ (eV) 


5.26 ±0.24 

9Q -i- n 98 
o.zy ie u.zo 


5.25 ±0.24 


5.27 ±0.24 

K 9E _|_ n 97 
O.oO IE U.Z I 




31.9° ±9.0° 

ZO.O IE i.O 


43.1° ± 15.1° 


50.4° ±23.3° 

7D 4- 97 9° 

i u.o ± Z( .z 


\C'\ (eV) 


3.71 ±0.50 

O.UZ IE U.04 


4.03 ±0.58 


4.15 ±0.51 

/l qc _|_ n AA 
4.00 IE U.44 


5 C ' 


-27.7° ±11.1° 

— lo.o IE ZU.O 


-21.4° ± 16.5° 


-17.1° ± 18.9° 

9 Q(30 _)_ 99 70 
— Z.OO IE ZZ. 1 


\P'\ (eV) 


40.0 ±2.6 
4y.o ie u.y 


32.6 ± 11.9 


30.6 ± 11.0 

9fi 9 4-8 9. 
ZO.Z IE o.O 




177° ± 1° 

7 no _i_ a oo 
I .y x 4.8 


9.04° ±6.95° 


-1.59° ±4.60° 

1 7Qo 1 pro 

— i/y ie 


-Pew,np ( e V) 


62.2 ±6.9 

iy.4 IE o.u 


21.1 ±4.7 




0EW 


-75.9° ±3.8° 
— yu.y ie o.u 


-87.2° ±6.0° 


- 


^EW,NP ( e ^) 


65.4 ±7.4 

9 Q« _L A OA 
O.yO IE 4.o4 


- 


86.4 ± 19.5 
«i 7 _i_ 91 1 

— 01. 1 IE Zl.l 


v'ew 


-| n /in 1 O O 

106 ± 3 
31.6° ± 78.3° 




146 ± 23 
27.7° ± 15.4° 


^np (eV) 




16.9 ± 12.2 


56.7 ± 13.3 
25.3 ±5.0 


0P 




-167° ± 16° 


6.91° ±6.93° 



-81.5° ±23.6° 



Table 6: The fits for NP (IV, V, VI). See the text for details. 
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NP fit VII 


X mi n/ d0 J 


7 


It 1 '! 
1 1 | 


r 


|o | 


Oc 


1 p'I 
r 1 




CO 1 in 


r 1 r\ 

0.0 ± U.z 


44 ± li 


4.1 ± U.O 


00 _i_ 1 n 

— zo ± iy 


00 ± 12 


D 43(Q3%1 


73 + 8 


5 2 + 03 


73 ± 8 


3 3 + 08 


— 15 ± 21 


46 ± 12 


2.6 (45%) 


69 ±31 


5.2 ±0.3 


33 ±41 


3.6 ±2.5 


-20 ± 25 


41 ±41 




P' 

r EW,NP 




p'C 




P' 

r NP 


0P 


179 ±2 


2.1 ±3.8 


9.8 ±52 


76 ±21 


-93 ± 22 


15 ±8 


104 ±37 


0.6 ± 1.3 


65 ±45 


33 ±29 


65 ±44 


-140 ± 30 


48 ±24 


48 ±38 


-172 ±8 


19 ±9 


89 ±7 


2.7 ±3.9 


-180 ±340 


8.6 ±42 


5.7 ±38 



Table 7: The fit for "NP fit VII". 



We list just three of them in Table 7. Since the Xmin' s °f "NP fit VI" are quite low, 
the solutions with low-lying Xmm values look similar to those of "NP fit VI" (See the 1st 
solution). The 2nd solution shows that all the 3 types of NP can be sizable. The 3rd 
solution corresponds to the "NP fit I" which is the NP in the electroweak penguin sector 
obtained in [5]. 

4 Conclusions 

We performed x 2 fitting to check if the conventional parametrization in the SM describes 
well the current experimental data of B — > nn and B — > i\K decays. Contrary to the data 
used in Ref. [8], the current data disfavors this parametrization given in (6) and (13) at 
2-3 a level. 

We interpreted this difficulty in the SM as a manifestation of NP and investigated 
various NP solutions. When a single NP amplitude dominates, NP in the electroweak 
penguin sector is the most favorable solution in accord with [5]. When two or more NP 
amplitudes exist simultaneously, solutions other than in the electroweak penguin sector can 
also explain the deviation very well. 
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